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' NOMEI3GIATI3RB 

All s 3 rnibols are defined whenever they appear 
first in the text. 

a ' - plate length 

b - plate width 

bg - effective width of plate 

- centroidal axis distance of stiffened element from 
mid plane of plate 

Constants equal to ^ 

G-^jyb., respectively 
d^ - spacing of stiffeners 

D^-jDy - flexural rigidities of plate in x-and y-directions, 

'z 5 

equal to E^h /(l-v v ), and Eh /(1-vv), respectively 

^ y y ^ y 

5 

l> 2 :y “ torsional rigidity of plate = G^yhv^ 

E - Young’s modulus for plate material 

E ,E - Young’s modulus for plate material in x- and 

y 

y-directions , respectively 
E* - B /E 

y 

E(x,y) - stress function, defined on page 58 

- shear modulus of rigidity 

h - plate thickness 

h -• equivalent plate thickness 

h„ - equivalent plate thickness (given by 13:801-1975) 

- moment of inertia of stiffener about the 
neutral axis of plate 
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^xy 

IT 

lT^,Ny - 

^xy 
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s - 

u,v,w - 
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a ; — 

®x>®y - 

Txy 


moment of inertia of the full area of multiple- 
stiffened element, including intermediate 
stiffeners, about its own centroidal axis 
torsional constant 

critical buckling stress coefficient 
numbers of half waves in which, plate buckles 
in X- and y-directions, respectively 
m7c/a 

bending moments per unit width 
twisting moment per unit width 
nix/b 

resultant normal forces in x-and y-directions, 
respectively 

resultant shearing force in xy-plane 
buckling load per unit width 
total compressive load in x-direction 
buckling or critical load 

perimeter of plate cross-section over width b 
displacements in x-,y-, and z-directions, 
respectively 

whole width of plate between webs (or from a web 
to edge stiffener) 

1/2 

arbitrary parameter = £ (P-Pq)/^^^^ 

strains in X- and y-directions, respectively 

shear strain in xy-plane 





A 

A 

O 

V 

V V 
x» y 

Cfy 


^xa 

c^max 

^cr 


- plate shortening in x -dir action 

- plate shortening in x-direction at buckling load 

- Poisson’s ratio 

- Poisson’s ratios for stiffened plate, in x-and 
y-directions, respectively 

- resulting stresses in x- and y- directions, 
respectively 

- resulting average stress in x-direction 

- maximum resulting stress along loaded edges 

- resulting stress at buckling load, in x-direction» 
_ mb/ a 



ABSTRACT 


The formula (13:801-1975) for the equivalent thickness 
of a multiple stiffened compression element, used in the 
design of cold— formed structural members, is reviewed. By 
carrying out an analysis for the buckling strength of an 
orthotropic plate, a more rational formula is suggested 
for design, 

Postbuckling analysis of a simply supported orthotropic 
plate is done using the method of successive approximation. 

In doing so, von Karman’s large deflection equations are 
converted into an infinite set of linear differential equation:, 
by expanding the displacement into a power series in terms of 
an arbitrary parameter. The curves of critical buckling load, 
load-shortening and effective width are presented for plates 
having different elastic properties. The possibility of a 
change in the buckling pattern of the plate is also indicated, 

Postbuckling analysis is also carried out by an 
approximate method similar to the one used by, for example, 
Ohajes (1974) for isotropic plates. 



CHAPTER I 


IHTRODUaTIOF 

1.1 Buckli.'ng. Postbuckling and Effective Width Concepts : 

Thin plate elements used in structures, are often 
subjected to normal and shearing forces acting in the 
plane of the plate. If these in~plane forces are suffi- 
ciently small, the equilibrium is stable and the resulting 
deformations are characterized by the absence of lateral 
displacements (w=0, u^^O, v O). As the magnitude of 
these in-plane forces increases, at a certain load inten- 
sity, a marked change in the character of the deformation 
pattern takes place. That is, simultaneously with the 
in-plane deformations, lateral displacements are introduced. 

In this condition, the originally stable equilibrium 
becomes -unstable and the plate is said to have buckled. 

The load producing this condition is called the critical 
load. The importance of the critical load is the initiation 
of a deflection pattern, which, if the load is further' increased 
rapidly leads to large lateral deflections and eventually 
to complete failure of the plate. - This is a dangerous 
condition which must be avoided. 

In the mathematical formulation of elastic stability 
problems of plates, one uses the neutral equilibrium assuming 



2 


a bifurcation of the deformations , That is, at critical 
load, of the possible two paths of deformation (one 
associated with the stable equilibrium, and other one 
pertinent to the unstable equilibria condition), the 
plate always takes the buckled form. In addition to the 
existence of this bifurcation of equilibrium paths, elastic 
stability analysis of plates assumes the validity of 
Hooke ' 3 law, 

Classic al buckling problems of plates can be 
formula.ted using (i) equilibrium approach which leads to 
an eigenvalue problem, (ii) energy methods such as Ritz 
method, Rayleigh's method, G-alerkin's method, (iii) numeri- 
cal methods using f inite-difference and finite element 
techniques, and (iv) dynamic approach. Methods (i)-(iii) 
have been utilized in solving buckling problems for 
isotropic, stiffened and orthotropic plates by many 
investigators; a comprehensive account is available in, 
for example, Ghajes ( 1974 ) and Szilard ( 1974 ). Tbe use 
of dynamic approach (Leipholz, 1970) becomes a must when 
the plate is subjected to nonconservative forces. 

Besides this classical buckling theory, the 
behaviour of flat plates after buckling is of considera- 
ble practical interest. The postbuckling analysis of 
plates is usually difficiilt, since it is basically a 
nonlinear problem and is markedly different from that of 
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thin struts. Vsfhile a small increase in the critical load 
for struts will produce a complete collapse, the critical 
load may be but a small fraction of the load that causes 
failure in plates. This increased load-carrying capacity, 
beyond the critical load, originates from the effect of 
large deflections and from the fact that the longitudinal 
edges of the plate are usually constrained to remain 
straight. The use of this additional strength is of great 
practical importance in the design of structures, since by 
considering the postbuckling behaviour of plates, considera- 
ble weight savings can be achieved. 

To study the postbuclcling behaviour of thin plates, 
one has to go for the solutions of von Karman’s large- 
deflection equations. Since the exact solutions of von 
Karman’s equations are cumbersome, generally approximate 
methods have been employed, for example, by Timoshenko 
and G-ere (1961) and Marguerre and Volmir (See, e.g.,0ha3es, 
1974 ). Wang ( 1949 ) and Brown and Harvey (1969) have used 
finite difference technique to solve large -deflection 
equations for plate subjected to uniform lateral pressure 
and compressive edge loading. Similar investigation has 
been made by Hooke (1969) using the perturbation method. 

The perturbation technique has also been applied by Stein 
( 1959 ) iu studying post buckling behaviour of isotropic 
plate subjected to longitudinal compression, Rushton(l970) 
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lias studied the same problem using the method of dynamic 
relaxation. For orthotropic plates, the only solution 
available is due to Prabhakara and Ohia (1973) who have 
solved von Karman's equations by expressing displacements 
and the stress function in double Fourier series. 


Till the buckling load is reached, the stress 
distribution along the loaded edges is uniform and' along 
the unloaded edges its intensity is zero. As the load 
is increased beyond the critical value, the distribution 
of stress along the edges becomes progressively nonlinear. 
Tests have shown that the loaded edges are more heavily 
stressed in the vicinity of the unloaded edges and that 
stresses remain virtually unchanged near the centre of 
the. loaded edges (Fig, 1-1 (a)). The most common method 
of accounting for nonunifom changing' stress distribution 
is to use the method of effective width, originally propo- 
sed, by von Karman (1932). 'When failure of the plate is 
impending, almost the total compressive force is carried 
by two strips, each of width bg/S, adoecent to the un- 
loaded edges (Fig, 1 -1 (b) ) . The quantity bg, the effective 
width is defined as 

1 b 



/ a .dy 

^max 0 


( 1 . 1 ) 


This implies that the centre portion of the plate is 
ineffective for the purpose of calculating the strength 
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of the plate. With this concept, the problem becomes 
one of finding the effective width b^ at a specified edge 
stress ^■''^hen bg is loiown, the following relation- 

ship s apply 


(i) 

(ii) 

(iii) 


o 

max 

®max 

‘^max 


m 

# 

CJ 

— 0 ♦ bQ=b , 

P= Or.Tr ^h 

cr 

^max 

^av * e ’ 

"^av 

cr* 

^max 


= %ax *>0*^ 

• 

J * 

^av 

^ult’ ^e~^ult 

’ “y '’ult*’ 


P T j. = ^ Tj-bh = o h T . h 

ult ult y ult 


In these I’elations is the average stress due to the 
applied total load P, is the value of at 

collapse, Oj is the yield stress of plate material, 

P■^3_^ is the value of P at collapse and. h^-^ is the value 
of bg at collapse. 


The effective width clauses in current Americal 
Specification (AISI, I 968 )* and the earlier Canadian 
Standard (GSA-S1 36, 1965) were based primarily on experi- 
mental investigations by Yfinter (1947), whose effective 
width formula is incorporated in these specifications. 
The current Canadian Standard (CSA-SI 36 , 1974) contains 
Q-n effective width formula of the type proposed by von 
Karman in 1932 . This change in the Canadian code was 



* IS: 801-1 975 specifications are also similar. 
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based on the work of Lind et al,(l976) who did a 
statistical analysis of the available experimental data. 

An excellent review of the experimental work concerning 
the effective width concept is given by Roorda and 
Y e nkatar amai ah ( 1 9 7 8 ) , 

1,2 Present Y/ork: 

Except the work of Prabhakara and Ghia (1973) > 
no theoretical investigation is available which deals 
with the postbuckling analysis of orthotropic plates. 
Moreover, the said work neither gives any relation for 
the effective width nor takes into account the effect 
of changes in the buckling pattern. In this thesis the 
postbuckling analysis of orthotropic rectangular plates 
of symmetric cross-section is carried out. Yon Earman's 
large-deflection equations, as applied to orthotropic 
plates, are reduced to an infinite set of linear differen- 
tial equations by expanding the displacements into a power 
series in terms of an arbitrary parameter. The first three 
equations of the infinite set correspond to small deflection 
theory whose solution yields the buckling load, vSolution 
of the suceeding equations gives the postbuckling solutions. 
The curves of loal-shortening and effective width are 
presented for plates having different elastic properties. 

The results obtained are compared with those of Stein 
( 1959 ) for isotropic plates* Postbuckling analysis is also 
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carried out by an approximate method, similar to that 
given in Chajes (1974). 

13:801-1975 formula concerning the equivalent 
plate thickness, for design of multiple -stiffened light 
gauge sections, is also reviewed, A more rational 
formula is suggested on the basis of the buckling 
strength of orthotropic plates. 



GH4PTER II 


MUITIPLB STIBBEIIED LIGHT GAUGE GOIffRESSIOH I/IEMBER3 

2 , 1 Introduction; 

Multiple stiffened compression elements are common 
in light gauge members made of steel or aluminium. The 
Indian Standard IS: 801-1975 specifications* for the design 
of cold-formed steel structured members specify that a 
flat compression element, to be considered multiple- 
stiffened, shall be stiffened between webs (or between a 
web and a stiffened edge) by intermediate stiffeners 
parallel to the direction of stress. If the stiffeners 
are spaced so closely that the sub-elements between 
stiffeners indiyidually do not exceed the specified limit 
of flat width to thickness ratio for fully effective 
stiffened compression elements, 'the flat width ratio' 
for the entire element may be computed as that of an 
’equivalent' fictitious, isotropic plate element without 
intermediate stiffeners whose width Wg is the whole width 
between webs (or from a web to the edge stiffener) and 
whose thickness hg is determined from 



* The AISI and OSA standard specifications are also similar. 



where Ig is the moment of inertia of the full area of the 
multiple-stiffened element, including the intermediate 
stiffeners, about its own centroidal axis. Infact the 
value of hg specified in Eq, (2.1) is the thickness of 
a solid plate element having the same amount of moment 
of inertia about its centroidal. axis as the actual, element. 

A closer examination of the mechanical behaviour 
of the element in compression reveals that the formula 
in Eq, (2.1) is based on a fallacy. Introduction of 
intermediate stiffeners transforms the plate into an 
orthotropic plate (Eig. 2.1) having much greater stiffness 
in the longitudinal direction than that in the transverse 
direction while the torsional stiffness is not very much 
different than that of the imstiffened plate. Since all . 
the three modes are present in the buckled element, it is 
clear that Eq. (2.1) overestimates the thickness of the 
equivalent plate. In this chapter, a more rational formula 
for the equivalent thickness of a multiple-stiffened 
compression element is suggested which is based on the 
analysis of the buckling strength of an ortho tropic plate, 

2.2 Basic Equations: 

To determine the critical in-plane loading of a 
flat plate, it is necessary to have the equations of 
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eqailibrium for the plate in a slightly bent configuration. 
For a loading condition consisting of constant biaxial 
compression forces (N^. and 1^) and constant in-plane 
shears (J^y) , this equation is expressed in the following 
form in terms of moments My, M^y and the lateral 
displacement w(x,y) of the middle surface of the plate 
(cf., e.g. , Timoshenlco (I 9 SI)): 


S^M. 


X 


2d\ 


■XJ 


2 

d M. 


+ 


6x' 


3xdy 


9y 


2 


2! + F 


-2 
3 w 


^2 
3 w 


+ H. 


3^w 


y 


ay 


2 ^ 2l^zy 


0 


3x3y 

( 2 . 2 ) 


Equations that express the moments in terms of 
displacements will be obtained by relating the moments 
to the stresses, the stresses to strains, and finally 
the strains to the displacement. For a thin orthopropic 
elastic plate, these relations can be written as 




D 

X 


0 


3^w/ax^ 

IL, 

y 


X y 

h' 

0 

< 

1 

1 

3^w/3y^ - 

^y 


0 

0 

^xy 

i 

2 

-3 w/ax3y 


( 2 . 3 ) 


where 


3„h5 




12(1- v^Vy) 




B h' 

y 




. a 


xy 


G- h- 
xy 


( 2 . 4 ) 
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Here, the plate is of uniform thickness h. The quantities 
5 are, respectively, the flexural rigidity per 
unit width of plate, the Young's modulus and the foisson's 
ratio in the x-direction while the corresponding quantities 


in the y-direction are By, Ey and torsional 

stiffness and the shear modulus G-^^y is expressed as 




3 ,+S +2 V E 
•^x y X y 


Further due to Betti's reciprocal theorem 


in which v and B are the Poisson's ratio and the flexural 
rigidity for an isotropic plate* Ifow, the substitution 
of (2.3) and (2.5) into (2.2) gives, the desired 



2.5 Buckling Strength* 

We idealize the problem by assuming that the 
edge fl an ges of the member prevent out-of -plane displacement 
of the edges of the plate. It is also assumed that the 
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flanges offer negligible resistance against rotation and 
in-plane displacement, Then the boundary conditions may 
be written as 

0 w 

N =1^ =0, w = — ^ = 0; y=0,b; 0< x < a (2.7) 

9y 

The ends of the plate are assumed to be restrained against 

out-of -plane displacements but force to move in the plane, 

and loaded by a constant compressive force of intensity h 

per unit length of mid surface: 

9 w 

N = - h, 11= 0, w = — — = 0; X = 0, a; 0<y< b 

( 2 . 8 ) 

Let us assume that the solution to Eq. (2.6) (with ]!^=l^y=0) 
is of the form 

“ °° m-rtx n-jcy 

w = E 2 A sin sin (2,9) 

m=1 n=1 a b 

in which m and n are the number of half-waves that - the 
member buckles into in the x and y-directions, respectively. 
The assumed solution satisfies the boundary conditions (2,7) 
and (2,8), and there remains only the task of ensuring 
that it also satisfies the differential equation. 
Substitution of (2.9) into (2.6) eventually leads to an 
eigenvalue problem; for a nontrivial solution we must 


have 
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a 2 m-n- 4 V m-rt 2 nu 2 

2(vD+D,,) {— ) (— ) 


"a 
nn; a 

- “y(- 


xy- 


( 2 . 10 ) 


To determine the least value of N, the expression on the 
r.h.s, of (2.lo) should be minimize w.r.t. m and n. It 
is found that h is minimum when n = 1 and m takes the 
value 


a D 1 A 

Consequently, from (2,10) we obtain 

V = ^ + I'xyJ 


( 2 . 11 ) 


( 2 . 12 ) 


For an isotropic plate, since 

D £h^ 

D = -IZ =B=: ^ 

^ Z 12(1-V ) 

(2.12) yields the familiar result: 

if 


(2.12a) 


2.4 Equivalent Plate Thickness: 


From (2.12) and (2.12a) it is obvious that 

V2 a 


1 ■^xy , A 

or 2 V- B , D 1 


(2.13) 
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It follows that the orthotropic plate has the same buckling 
load as of an isotropic plate of the same material and 
width and having the thickness h given by 


h = h [ 


1 D D 1 /2 

1 ((-^) ' 

2 d2 


4“ V 4“ 


^xy ,V3 


(2.14) 


For a stiffened plate (Fig. 2.1), various rigidities 
can be evaluated in the following manner (cf . , e.g.,Ssilard 
(1974)): 


Eh' 


D 1 B = 

y 


B, 


EhCy 


1 2(1-v2) 
2 


El 


B + 


1 - 


+ 


n 


x; 


dv 


J. 


(vD+B^y) 1_B + 


X 


4(l + v) 


(2.15) 


In (2.15), is the torsional constant, I^^ is the moment 
of inertia of stiffener about the neutral axis of the plate 
and c^; j ix are explained in the Fig. (2-1 ). u.se of (2.,15) 
in (2.14) results 


h 

h 


1 




h 


iyh3 


) +1+3(l-v) 


Jx 


d2^.5 
(2.17) 

Alternatively, the expression for equivalent plate thickness 
can be derived by evaluating the rigidities in a manner it 
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is done for corrugated sheets. Denoting the perimeter 
of the plate cross-section over the width b by s (following 
the loops of the stiffeners), we, therefore, represent 
rigidities as 


b , I3 b 

Tj V — Pi "d D, D — ( 1 - ^ )D 

b b 

V D v — 'I'D, (v D +D ) = — D 


X y 


X y xy' 


(2.18) 


Use of (2.18) in (2.1 4.) and a subsequent simplification 
yields 


h 

h 



<15 i’ 


(2.19) 


Equations (2.17) and (2.19) are suggested as 
replacements for Eq. (2.1) to determine the 'equivalent' 
plate thickness in light gauge design. Numerical calculations 
show that there is a substantial difference between the 
value of h determined from(2.l) and (2.17) or (2,19); 

Eq, (2.1) really overestimates the thickness of the 
equivalent plate. 



CHAPTER III 


POSTBUOKLIHG Ax^!ALYSIS OP ORTHOTROPIG PIATES-I 

3 • 1 I ntroduction: 

A number of studies are available in literature 
concerning postbuckling analysis of flat rectangular 
isotropic elastic plates, tlie most comprehensive being 
that of Stein (1959). Except the work of Prabhakara and 
Ghia ( 1973 ), there is no investigation which deals with 
the post buckling of orthotropic plates under compression. 

In this chapter, the postbuckling analysis of rectangular 
orthotropic plates of symmetric x-section is carried out. 

A method of successive approximation, originally applied 
by Alexeev (1956), is utilized. Yon Ear man's large deflec- 
tion equations are converted into an infinite set of linear 
differential equations by expanding the displacements into 
a pov/er series in terms of an arbitrary parameter. The 
first few equations of the infinite set turn out to be 
the small deflection equations. Solutions of these and 
succeeding equations permits a study of the behaviour of 
the plate at buckling and beyond, up into the large defle- 
ction range. The curves of load-shortening and effective 
width are presented for various cases of plates having 
different elastic properties. The possibility of a chaaige 



in the buckling pattern of the plate is also indicated 


3 • 2 Q-overning Equations : 


For a plate with no lateral load, the von Karman 
large-deflection equations can be vn^itten in the form 


i^x,x 


= 0 

(3.1a) 


^y,x 

= 0 

(3.1b) 

Vw - 

(M,xx+^^ 

w +2^J w ) = 0 

T ,yy xy'",xy ^ 

(3.1c) 

where ^ w = 

= D w 

+ 2 ( V D+D . )w +D w 

^ ^^■^xy^'^jxxyy y ,yyyy 



Here a comma denotes partial differentiation with respect 
to the suffix that follows. 

'The strain-membrane force relations are 

(3.2a) 

(3.2b) 

y ' = _ (5.2c) 

h G-xy 

The forces in Eqs. (3.2) can be solved for and are 
expressed as 


= 


By - 


1 

(* 


:x: 


^x 


(i 

h By 


V 

y s 

j 

V — ) 

^ T? 

Jlrv 


1 N. 
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IT, 


E h 

V 


y 


1 -V V 

^ y 






The strain-displacement relations are 

1 o 


u - w 

2 > 


X 


1 2 

-y- + - W 

>y 2 »y 


(3.3b) 

(3.3c) 


(3.4a) 


(3.4b) 


'xy 


= U^^ + Vv+W^W 


,X ,y 


(3.4c) 


It is assumed that u, v, and w may be expanded in a power 
series in terms of an arbitrary parameter a: 

(3.5a) 
(3.5b) 
(3.5c) 


0,2 4 

u = au + auo+ a u,i+ 

o 4 

o 2 4 . 

Y = a V + a -Vq + cc V.+ 

o 4 


3 5 

w = aw-j + a + a w^ + 


Por plates (without initial eccentricities) subject to 
in-plane loading, the deflection w is zero in the loading 
range prior to buclcling but u and y haYe values other than 
zero. The series for u and v is therefore expected to 
start v\rith zero power of a while the series for w is expected 
to start with a nonzero power. Prom the available solutions 
of the postbuckling behaviour of rectangular plates subjected 
to longitudinal compression, it is seen that the square of 
the, centre deflection is nearly a linear function of the 
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applied load. For the compression problem it is convenient 
to relate a to the load P so that a = (P-P^)/Pq, 
where P^ is the critical load. Thus, the series for w 
will start with the first power. The odd powers in the 
series of u and v and the even powers in the series of w 
vanish for problems of the type considered and, for 
simplicity, they have been omitted from the start. 

It is advantageous to expand the externally applied 
loads also in terms of the arbitrary parameter' a. Upon 
substitution of (3.5) into (3f3) and (3.4), the following 
relations are obtained: 


■xmn 

m= 1 ,3 n=1 , 3 


n=0 , 2 


OO OO OO 

U = S u E E 

^ n= 0,2 ni=l ,3 n=1,3 

OO CX3 CO 

^xy ~ ^ "Scyn “ + ^ 2 ^^xymn°^ 

n=0,2 ^ m=1,3n=1,3 


where 


II 

*^1 ^^n,x '^y^n,y 

II 

^2^ ^x^n,x‘^ 

1 ■ = 
XJIi 

^3^^n,y'^''^n,x^ 


(3.6a) 

(3.6b) 

(3.6c) 

(3.6d) 
( 3 . 6e ) 
(3.6f) 
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0 


1 


(w yj + V 

xmn 2 ni,x 'n,x 


2 ^ V^m,x'^n,x 


^^xymn ^3 


w w ) = N 

y m,y xma 

(3.6g) 

+ w w ) = M 

m,y n^y-^ ynm 

(3.6h) 


(3.61) 


in which. 


Gl = ExVd-v^Vy) , 02 = EyJl/d-v^'-y), O3 = 

Since a is taken to be an arbitrary parameter, the 
stipulation that a power series in a vanish requires 
that each coefficient of power series vanish." Substituting 
(3.5c) and ( 3 . 6 ) in (3.'!), this requirement leads to the 
following linear equations, which are the first few of 
an infinite set: 


^xo,x ^yo,y 

w + W = 

yo>y -^xyo,x 

VW1 -(Exo*1,xx'^ V 


0 

0 


w 


I >yy' 


.+ 


) 


^■"^xyo ,xy^ ® 


(3.7a) 

(3.7b) 


^2,x ^xy2,y “ ^®x11,x ^y11 ,y^ i 

-V,y ^'^y2,-x = - ^"\ii,y 

Vw^ -(IS^Q ^3, XX ^yo^3,yy ^®xyo ^3,xy^ 
= (Hx2+i^x11 ,xx ^'l »yy 

^ ^^xy2 ^y1 1 ^ ,xy 


(3.7c) 


,.('^.7d) 
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^4,x'*'^y4,y 5,x'''^y1 3,y'^^y31 ,y^ 

®y4,y'^^xy4 jx"^" 3 ,y'^^xy1 3 ,x'''^xy3l ,x ^ 


(3.7e) 


^^5"^^xo^5,xx+V^5,yy''2Vo 

~ ^^2’*’^11^ ^3»xx ^^y2'*'^y'il^ '^3>yy 

^^y2'^^xy1 1 ^^3jXy 3^^’1 ,xx 

+ (Fy 4 + 2 I^^ ^)w^ ^yy ^ ^ ^ ^ y 4 ^^ ^ y 1 3 ^ ^ y 3 1 ^ 1 , X y 


(3.7f) 


3.3 Solution Procedure; 

llie problem to be solved is the post buckling behaviour 
of a simply supported plate under longitudinal compression 
with edges constrained so that the displacement of each edge 
in the plane of the plate is uniform. Accordingly, the 
following boundary conditions are imposed on the plate: 



Zero deflection: 

w(o,y) - w(a,y) = w(x,o) = w(x ,b) = 0 (3,8a) 
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Zero moment : 

w^x^(o,y) = w^^^(a,y) = w^yy(x,o) = w^yy(x,b) 
Constant displacement; 


= 0 

(3.8b) 


^jy(o,y) = u^y(a,y) = v^^(x,o) = v^^(x,b) = 0 
Zero shear stress: 

"^jX = ^,x(^>y) = u^y(x,o) = u^y(x,b)= 0 

Loaded edges: 
b 

i<-V^=o,o. 4y = 


(3.8c) 


(3.8d) 


(3.8e) 


Unloaded edges; 

a 

•^(U) „-h'i2: = 

o y y=o,t) 


0 


(3.8f) 


2 

Substituting (3.6a,b) and noting that P = Pq+oc P^ , the 
preceding force boundary conditions take the following 
form for different approximations in a : 


loaded 

b 

/ 

0 

b 

/ 

0 

b 

I 

o 


edges: 

^^’^xo^x=o,a/^^ ~ - ^o 

^^"^xZ ^Sc1 1 ^x=o,a “ " 
^^X4 ^^x 1 3 ^x-o ,a ^ 


(3.9a) 

(3.9b) 

(3.9e) 
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Unloaded edges: 

3» 


/ 

0 

^~yo^y^o,b ^ 

(3.10a) 

a 



f 

0 

(17^2+^^-! 1 )y^Q^b = 0 

(3.10b) 

a 



/ 

0 


(3.10c) 


3.3.'! Fundamental Solution: 

Eqn. (3.7a) can be written in terms of the 
displacements u^ and : 

t ^o,xy^^2^o,yy^^5^o,xx “ ® 

The solution of the above equations satisfying the boundary 
conditions (3.8c,d) and 3.9a) is 


u. 


fi bh 

X 


(x - I ) 


^ _ - ) 
o S^bh " 2 


It therefore follows that 


(3.11a) 


(3.11b) 




» ^yO ^'^yO ® 


^’xo 

Using (3.12), Bq. (3.7b) can now be written as 


(3.12) 
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Vw^ 


+ 



w 


1 ,xx 



(3.15) 


whose solution can be taken to be of the form 


mn; nTt 

Sin rix Sin . >Ty; M = — , . J = — — (3.14) 

a b 

v\?hich satisfies the boundary conditions (3.8a,b). 

This solution req_uires that 


Pq = bM"^[ D^M^+2( j (3.15) 

So far, the solutions obtained are identical to 
those obtained from the small -deflection theory; Pq, 
therefore is to be identified as the buckling load. The 
lowest value of the buckling load is determined by the 
choice of rn and n for a particular aspect ratio a/b. 

It may be noted that the amplitude A-| can not as yet be 
determined. 


3.3.2 First An-proximation: 

Sqs. (3.7c) are rewritten in terms of U2 j''^ 2’ : 

'’y°U°3^'^2,xy = " .xx™! ,x'*'''y"l ,xy’'l ,y 




TS" 


(v^a2+05)u2^^y+a2V2,yy+<33y2,^x 

= -°2f"i,yy''i,y'" ,xy'"i ,x > 

.xy™! .x+^l .y"! ,xx) 
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Upon substitution for w-| from (3.'i4) into preceding equations, 
the following solutions for up and vp are obtained. 


Up = - r + — 

^ ^ E^bh 8 


A? o 


, 2 ,2 

p a A^ M -V U 

Ml (x- -• )- ~ r ^ Sin 2Mx 

■12 16*- M 


-M Sin 2m Oos 2 U 7 ] 


(3.16a) 


■^2 = [ 


P V 
O X 


9 b A -1 

-n ( y - 2 )- ^ [ 


s^bh 8 

-K 


U Cos 2m Sin 2% ] 




Sin 2Uy 


(3.16b) 


So that 


^ x 2 + ^\11 


V ' V11 


P. E,^h 2 

— - -2^- — - M Cos 2 % 
b 8 

E h p 

.JE — # Cos 2MX , 


(3.17) 


V 2 -" V 11 = ° 


how must be determined from Eq. (3.7d), After 
substitution of the IT's and w-| , ( 3 . 7 d) becomes 

p 2 

Vw. + ~ w_ = A.[-^ M _ ~ h(E M^+B iV^)] Sin m Sin By 
3 b 3, XX i>-^ Ig ,x y J 


Eh . Eh Af 

+ 1 Sin Mx Sin 3% + 

1 6 S 


Sin 31'is. Sin Ny 
(3.18) 


It should be noted that Sin Jfe. Sin Ey is a complementary 
solution of Eq. (3.18), and a term of this kind appears 
on the right hand side of this equation. No solution to 
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-Sq. ( 3 . 1 8} is possible that satisfies boundary conditions 
of this problem unless the coefficient of this term on the 
ri<ght hand side of this equation is zero. Thus, 

p 

This relationship provides the means of determining the 
value of . If the trivial case = 0 is ignored, the 
amplitude A-| is found to be 

^2 ^ ^ 4 ^ s (3.19) 

' bh X j 

Now the complete solution for which satisfies ( 3 . 18) 
and the boundar^r conditions nan be written as: 


W 2 =A^Sin Mx Sin Ny + A-| ^Sin Ik Sin 3%+-^5-l Sin 3i&xS±n Ny 

(3.20) 

where A^ can not be determined as yet, and 

^ M^[ D l/+2( vD+D )M^(3N)^+I> (3N)^- — ]“^ 

‘2-15 ^ y b 

h A^ I* 

A., = ^ ^1^ [D^(3M)'^+2(vB+D^y)(3M)^N^+D^^N^- ^(3M)^] 

1 6 


2-1 -1 


' 3 I 


y 


3 . 3.3 Second Approximation: 


I'irst the expressions for and 's can be 

written in terms of w^. Then, Eqs. (3.7e) become 
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“l"4,M'"°3"4,yy+(Vl^3K,=cy = 

-(^“1-^3) ("'1 ,xy"'3,y"'*1 ,y'*'3,xy)-03("l ,x™3,yy+*3,x'"l ,yy) 

( ''x02+Gj)u^^^y+03y^_^^+C274.,yy 

= - °2('"l,yy'"3,y+'''l,y’»3,yy)-{ Vg+Cj) . 

(”1 ,xy’"3,x+'"l ,x'’3,xy)- °3("1 ,y''3,xx-"’'3,y’'l ,xx 


) 


The solutions that satisfy the boundary conditions are 
2 

•1 "1 ^ Cl ’A. "lop P P 


U4 = 


A A 

Sin 2m - V Sin 4Mx 

1 6m y 

P (Gpl^^-v 0. I:^) (A-j ,+A:,h ) 

[ A_-A4 ,+3A,-, -40,ir — — -r-^- — J — 

8 ^ 5 13 31 3 0^051^+020^11'^+ err 

A^A.,M p (402^^- V G.M^) 

Sin 2l\feCos 2^j + ^ [I+O^r ^ 


+ 


A^M 


8 " ^ 0 ^ 0 ^M^+O 2O 5 ( 2 F) '*•+^1^ ( 2 IT) 

44jy--|- U [1 _ O5M . 

( 40-, M^+0-,02(^~^xV^^/^3-^x^2^^ 


,] 


Av|A^>iM p 

Sin 2m Gos 417+ • — — [1 - G^’' 


.] Sin 4 Mx Gos 21^ 
0^05(21)^0205^0(21)^1^ (5 2^a) 


where 0 = 0^20- ~ ^ l'^3''j~^2'~'3^ x 


A^A5K^ 


4 

-] A 

iTn 


v..v_P - G.G^v^^-OpO^v. 

,2 . , ,2 


-__liL(v- - )- h rA.(U^-v M^)+A^,(3/+\Mt ] Sin 2% 
T4 — 2 ^ 8U L 3 X 13 X J 


A*^ A-| ^ 


2 

- )Sin4lfy+ “^[ "^3-^31 ■**5-^1 3 


P (G^I — V GpH ) (A-i •v'^'A't-i ) 

- 4'G3r — y . JA 3, j Gos2Mx Sin 2% 

04051^+02051^+^1 IT 



50 


+ -1— U — [l-O^N^ — ;! — ] — I-_l_ '] Go s2Bf!&: Singly 


+ -J— 1L.[1+C3I 
8 


( 4021 ^^+ 0 -, G 2 ( 1 -V ^)M^/G,-v ). 

f ^ J.M l■.l . l.■ - ^y ^. . .I . , , . .. . - jiZ 

^ O^G^M^+G 203 ( 211 )^+ 0142 ( 21 ) 2 

(4G^m2_ ^ q^jj2^ 

^,^2 — ^ — 1 Cos4I&Sin2l'Ty 

2 G^G^(2M)'^+C2G^ir+G(2M) " 


ir 


(3.2lb) 


So that 


A-,M 


N^^+21jj^3= ^ [ (A^-A^ 3 )Gos 2 Ny+A-| 3 G 0 s 4]5ry 


+4G2G^i^(A^^+A^^ ) { G^G^MVG2a^lJ^+GM^.^ r”^ Go s2ihJ:Gos2I’ly 
_4C2G3iV^A^ 3{G^G^M^+G2G3(21 J)Vgi/( 21T)^ r’'oos2M:!i:Oos4Ny 
_G2C3A3-,{G^a^(2M)^+G2G^l^^+G(2M)^lI^ Co s^i^xCo s2l^yj 




li 

^ -Y- 


( 3 , 22 a) 


[(A^-A^^ )Gos 2 Mx+A^^Cos 4 ivK 


+40^0 (A^ ^+A3^ ) ^ G ^ C ^M^+G20^1T^+GM^1^ Go s214x Go s 2 Ny 

_G^Gyl\^^{G^G^xM'^+G 203 ( 2 N)‘^+Cl/( 2 lO^ 0os2aJxGos4% 
4^ {G G (2^^'^^''^ ^ /■^T^4',^/on/^^^ t\t2 1"“^ 

3.t3 


-40^0 3M'^A^:, { G ^ G ^ (2M) \G2G3 ( N) ^+C ( 2M) }“ Go Oo s2 ITy j 

(3.22b) 


A . M^U- 

N .+N ,^+l ,^= - G„(1- V V ) r 2G.Go(A.,+A,i ) 

^^xy 4 xyl3 xy 3 l 2 5 2 c y'^ l 12^13 31 ^ 


■ { 0^0^M^+G^G^+0¥^'^r^ Sin21CxSin 2Hy 
- G^02A^3 {G^G^M^+G 203(2^)^+014^(21)^}“^ Sin2l\!x Sin4l|r 
{0^GAm)^+0^.^+^ f'^ Sin4l4x Sin2lfy j 


(3« 22c) 
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1'iie differential Eq. (3.7f) for is now completely 

determined except for A^. The constant is determined in 

the same way in which A^ was determined: thus 

1 6P -1 

(3.23) 




3bh Aljf 

Therefore, upto second approximateion, the displacement 
functions are 


0 2 4 

u = a Uq + a u^ + cc u^ 

0 2 4 

V = a + a ^2 + a 


w = aw-^ + aw. 


(3.24) 


With a = (P_P^)/Pq 


3.4 Main Results : 


Some of the important results can now he written 


down: 


Oritical Load ; Using the expression in (2,4) for various 
rigidities, the formula (3.15) for critical load can be 
written as 




^o = 


. K 


i Where 

^ ^ mb 2 

K = (-) +2V 


4(l-v^Vy) 


a ^2 

y E,/E +1+2V, \ 


' X 

is the nondimensionalized buckling coefficient, 


(3.25 ) 
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Load-short eniioff: relationship: Total axial shortening A 
is the sntn of the inward displacements at each end. Since 
u is positive in the positive x-direction» 


A = u(o,y) - u(a,y) (3.26) 

Substituting the expressions for nQ,U2, and u^ in the 
first of (3.24) and then using (3.26) will yield the 
following relation after lot of simplifica,tions: 




P P ^ 2 p 2 r 5('I+2Ki+P) 

/ + ( 1 ) +( 1 ) [ — — 

^ Pq ^o 'I+P ^o 4(1 +P)^ 

.2 


(— ! + — )}] 

2(K^+5p) 9-P 


(3.27) 


where is the axial shortening corresponding to critical 
load, i.e., Aq = ^Qa/B^bh 

and p = 




mb 

X 


K, 


V + 

y 


2(l_v V ) 
' ^ JJ 

1+E^A+2v^ 


'X/ y 

Bff ective-width relation: The effective v^idth can 

conveniently be defined as ■ 

P 


_A_ 

a 




(3.28) 


Substituting Sq. (3.26) in the above expression, we get 
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b 2 F P P P 

[_ +2( — _1)/(i4-^)+(^ _1) {3(i+2K^+p)/ 

D o ^ ^0 

4(H-p)^{1/2(K^+5fB)+pV(9-P)}}r^ (3.29) 

3.5 Discussions : 

Numerical calculations for the results are performed 

for vaocious plates having different elastic properties 

defined by the parameters E„,iil , v v . However, for the 

z y A y 

purpose of presentation here, the following three cases 
are considered; 

(i) S,^/Ey(tb be denoted, henceforth, by E*) = 1 , 

= Vy = V = 0.3. It corresponds to an isotropic 
plate. 

(ii) S = 4 , = 0,3 , Vy = 0.075 

(iii) = 0.5 , v^= 0.15, = 0.3 . 

In Bigs, 3-1 to 3-3, the critical stress coefficient K 
is plotted for three cases mentioned above. The portions 
of the curves defining the lowest critical values of the 
load are shown by full lines. The different a/b ratios at 
which the transition from m to (m+1 ) half-waves takes place 
for a particular value of E can be read of from the 
corresponding curve. Bor a given aspect ratio, the half- 
wave at the instant of buckling are more for B < 1, as 



compared with the isotropic case, and are less for 3* > 1 . 
For example, an isotropic plate of a/b ratio 3 will buckle 
in 3 half -waves whereas plates with B = 4 and 0.5 of the 
same ratio v>/ill buckle in 2 and 4 half-waves, respectively. 
Therefore, the value of the parameter ^(=mb/a) in these 
cases are 1 , 2/3, 4 / 3 , In Fig. 3 - 4 , the variation of K 
is plotted for various values of X ; it monotonically 
increases with,., X for E* > whereas for E < 1 it first 
decreases and then again increases as X is increased. 

Figs. 3-5 to 3-7 represent load shortening curves 

for a/b = 1,0, 1.5 , 5.0 respectively. For each a/b ratio, 

* 

all the three value of B are taken, load-shortening curves 
are drawn using the relation ( 3 . 2 ?) in which the values of 
m correspond to the lowest bucklir^ load; n always equals 
unity for this proble.m. In addition, the curves are also 
given for other values of m which intersect with these 
basic curves for the range plotted. There are similar to 
those obtained by Stein (1959) for isotropic plates. The 
intersection of the load-shortening curves indicate possible 
changes in buckling pattern. For finite plates, changes in 
buckling pattern are often observed experimentally. However 
experiment by Stein (1959) shows that the change does not 
necessarily occur at the intersection points of these lines 
for different numbers of halfwave,. 
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Load-shortening curyes also give an idea of the 
postbackling strength of the plate which depends upon the 
ratio 3 and the parameter X . Fig. 3-8 (a-c) show such 
curves for^= 1.0, 1.5, 2.0 respectively. For X = 1, all 

the plates show distinct postbuckling strength; plates 

* 

with ill <1 have a higher strength compared to plates 

with 3 >1. Ihe situation is .just the reverse for X= 2, 

For 1.5 anh 2, the postbuckling strength of plate 

with S = 0.5 is not much different than that of the 
isotropic plate. 

Fig. 3 “9 is a plot of non-dimensional effective 
width (bg/b) versus ^ for PA'q = 2, based on .iiq.(3.29). 
Figs. 3 -10 (a-c) , on the other hand, give the variation of 
bo/b with nondimensional load P/P^ for X=1.0, 1.5 , 2.0, 
respectively. In Fig. 3-'^0(a) the dotted cuerves correspond 
to the results based on the first approximation only. If 
a/b is given, the designer can determine the appropriate 
value of m from load -shortening curves, for a selected 
value of P/pQ and the known properties of the plate. This 
determines X , which is mb/a. Once X is known, the 
corresponding value of the effective width can be found 
out from Fig. 3-9 or 3-10. 



F 19.3-1 Buckling stress coefficient H (En I Ey»1 ) 
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A^Ao 

Fig.3-5(b)Nondimensionat load short«ning curv« 

(a/b =1 & Ex/Ey*4) 



Pig. 3-5(0) Nondimenstonal load shortening curve 
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OHAPTiJR IV 


POSTBUGKLIVG AIIiLYSIS OF ORTHOIROPIC PIATES-II 


4-.'! Introduction; 


A fairljr rigorous analysis of the postbuckling 
strength of orthotropie plates has been carried out in 
Chapter 3. The aim in this chapter is to present an 
approximate analysis of the same problem. Results are 
somewhat less accurate to the extent that they match with 
the results obtained from the 'first approximation' only 
in Chapter 3. The method consists in assuming a simple 
function for w (the mid surface displacement), solving 
for w and a stress function F(x,y) and then evaluating 
the arbitraxy constant in w by means of the Gal er kin's 
method. A similar approach for isotropic plate can be 
seen in Ghajes (1974). 


4* 2 Governing Relations: 

As stated in Chapter 3> the equations of equilibrium 


are 


811 

8x 



9^y 

9y 



3X 


0 

0 




(4.1a) 
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9 w 




w 


9x 


+ D. 


.2 

5 w 


■N. 


y ey2 


- P 




a2 

9 w 


9x9y 


d^w 


9y4 


0 


N. 


X 


0 w 
Sx 2 


(4-. 1b) 


The middle -surf ace strains are expressed in terms of 
the middle -surf ace forces, in the following manner: 


mmm 

( . 

h 

Ex 

1 



- ( 

h 

E 

y 

1 

%y 


V ) 


y B 


J 


X 


A 

E„ 


h *^xy 


/ 




(4.2) 


Three equations in (4.1a,b) involve four unknown functions 
in x,y. Hence a fourth equation is needed to determine all 
the functions. It is obtained by making use of the compa- 


tibility equation; it is a relation amongst ^niddle— surf ace 
strains and displacements and is of the form: 


2 2 

a 3 

£ + 2. 

0y2 ax^ 


2 

a r 


EL 


6x 9y 


3 w 2 

5x9y 


.2 

9 w 3 w 

3x^ ^ 


(4.3) 


To reduce the number of equations that must be solved 
simultaneously, a stress function l'(x,y) is introduced. 
Eqs, (4.1a) will be satisfied if the in-plane forces are' 
defined, as under 
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2 

6 F 

\ ^ TT 

9y 




’ ° 3^2 ’ 




h- 


9x-3y 


( 4 . 4 ) 


Substituting the relations (4»2) and. (4.4) into Eq. ( 4 . 3 ) 
and the relation (4.4) into Eq. (4.1b), we obtain 

6^ 1 0% 


1 a'^'p _ 2 _ 2 v^ ^ ^ ^ 

2 2 
0 w 0 w 


=(■ 


dx^ ay^ 


0^w 2 

Qxdy 

( 4 . 5 ) 


and 


3t 


+ 


h 0x4 

„2„ .2 


dj'^ 9x' 


2(vD+Dxy) 

^ 4j_ 

a w 


a'^w 

h 

2 . 2 

dx Oy 

h 

‘ 

2 

^ a p 

.2 
a w 

4 - 

a^p 

^2 

! W 

■ ax 3y 

8x3y 

ax'^ 

ay2 


( 4 . 6 ) 


We have thus reduced the number of equations 
that must be solved to two. Sqs. (4.5) and ( 4 . 6 ) are 
in the modified form, the von harms-n large deflection 
plate equations. 


^ • 5 Approximate Solutions: 

Let us consider the simply supported rectangular 
plate subjected to a uniaxial compressive force Since 

a finite-deflection analysis involves the deformations in 
the plane of middle surface as well as transverse bending 
deflections, both in-plane and transverse boundary conditions 
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must be spacified. The trajisverse boundary conditions 
corresponding to simply supported edges, i.e. 


.2 

5 w 

dx2 


x=0 ,a 


a^w 

W = — r = 0 
ay2 


y=0,b 


(4.7) 


With regard to the in-plane boundary conditions, the 
following assumptions are made: 

1 . All edges remain straight and the plate retains 

its rectangular outline during bending. 

2. The shear force vanish along the four edges 
of the plate. 

3. The edges y = 0, a are free to move in 
y direction. 

These boundary conditions are similar to ones stated 
in Chapter 3. 

In the present case the edges are permitted to 
move, provided that they remain straight. Thus only the 
average of Ny, and not Jly itself, must be equal to zero, i.e. 
a 

I dx = 0 

similarly (4.8) 

I dy = - Jfb = - 

Where is the average stress intensity along x = o,a and 

the negative sign indicates that the applied load is 
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compressive. 

The first step in the analysis is the choosing 
of a suitable function for the lateral deflection; we 
assume that 

mitx nity 

w = f Sin Sin — ( 4 . 9 ) 

a b 


Having assumed the form of w, we next solve Eq, (4.6) 
for the stress function E, Substitution of ( 4 . 9 ) into 
(4.6) gives 


1 

By dx^ 


1 

G-xy 



1 a4y 


f mrt 2 n-n; 2 Tmixx 

5 -{-) (- ) [Cos — 


0o)S 


2 n-n:y 


(4.10) 


The solution of this equation consists of complementary 
and a particular part. That is, 

j; “C p 

To obtain the complementary solution, one sets the r.h.s. 
of the equation equal to zero. But this is equivalent to 
letting w = 0 , Thus the complementar^r solution of (4.I0) 
corresponds to the in-plane stress distribution that exists 
in the plate just prior to buckling. At that instant the 
in-plane stresses are known to consist of a uniform stress 
F^(=N) and Hy=l'Cxy =0 . 


Hence 
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= Ay^ (4.11) 

In view of the first relation in (4.4) and the fact that 
~ ^xa^» can be rewritten as 



( 4 . 1 1 a) 


Ihe particular solution can be of the form 

2mmx 2mty 

= 3 Oos — + 0 Gos -r — (4.12) 

Pa b 


Substituting this expression into ( 4 , 10 ) and equating 
coefficients of the like terms, we obtain 


B = 


!£ 

32 


an 2 

(— ) and 
^bm 


e = -2- 


32 


bm 2 


Thus 


P. 


P 32 


• stn 2 2m7tx ' 

.(— ) Cos + E. 


y mb' 


X 


bm 2 2 n 7 ty-, 

(-_ ) Cos 
an ^ b 


such that the complete solution of 3q. ( 4 .I 0 ) is 


f , ,an 2 „ 2mTix ^ bm . 2 2mxy 
32 y mb^ a X an b 


"^xay 


(4.13) 


In order to find out the value of f in Sq. ( 4 .I 3 ), which 
should also satisfy Eq. ( 4 . 5 ), it is necessary to solve 
Bq. ( 4 . 6 ). Solution of ( 4 . 6 ) by the Galerkin's method 
gives us (cf.,e.g., Chajes (1974)): 
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^ ^ Effective Width Relation; 

As the buckling load is increased beyond the critical 
load, the distribution of the stresses along the edges 
becomes progressively nonlinear (see Chapter 1). The loaded 
edges are more heavily stressed in the vicinity of the 
unloaded edges and the stresses remain virtually unchanged 
at the centre of the loaded edge, where the stresses are 
equal to (or less than) the critical stress 
effective width b^ is, therefore, defined as 


b 


e* 


a 

max 


b 

/ a dy 
0 ^ 


Substitution for from (4*10), the above relation 
yields? 


E, , ^cr , ^cr 


X_ (1 


5B E 

^ X ^y 


o. 


) + 


max 


a ■ 
max 


(4.19) 


Or, in terms of the critical load Pq and the applied load P, 
b^ 




(3E^^'^+Ey)(l-PQ/P)+P^/P(E^XV Sy) 


(4.20) 


The form of Eq. ( 4 . 20) is similar to the one derived in 
Chapter 3 Eq. (3.29) . "When the ratio bg/b is computed 

for a given X (mb /a) and the ratio Bx/Ey, the values 
exactly match with those obtained using 'first approximation’ 
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in Gliapter 3; these values are shown by dotted curves 
in Fig. 3-l0(a) • Hence, to a certain degree of approxi- 
mation, the method employed in this chapter is similar to 
the method of successive approach presented in Chapter 3. 



CHAPTER V 


COHGLUEIHG RBIVIAPJCS 

This thesis has been concerned with the analysis 
of buckling and postbuckling behaYiour of orthotropic 
rectangular plates under longitudinal compression. First, 
a mult iple -stiffened light gauge compression member has 
been idealized as an orthotropic plate, simply supported 
along the edges. By carrying out an analysis for the 
buckling strength, a new formula for calculating the 
'equivalent plate thickness* has been suggested for design. 

It turns out that the 13:801-1975 formula grossly over- 
estimates the equivalent plate thickness. 

Y/hether it is a mult iple -stiffened light gauge 
section or a fibre reinforced composite plate, each can 
be treated as an orthotropic plate for the purpose of 
analysis. Postbuckling analysis of such plates has been 
carried out, first by using the method of successive 
approximation v/herein the displacements are expanded into 
a power series in terms of an arbitrary parameter, and 
secondly by an approximate method wherein the differential 
equation is solved using the G-alerkin's method. For these 
plates, an effective width equation has been derived and 
curves for load-shortening have been obtained. These curves 
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show a possibility of a change in the buchling pattern of 
the plate. 

A linear set of equations has been deriyed to 
replace the nonlinear large-deflection equations for 
plate and is shown to have the advantage of simplicity 
of solution, since much more is known about solving 
linear partial differential equations than about solving 
nonlinear partial differential equations. However, the 
linear set of equations are subject to certain limitations 
depending upon the application desired. It is to be 
expected that solution to certain problems might not 
converge satisfactorily and it appears that the linear 
equations cannot be used to solve postbuckling problems 
for plates with initial eccentricities. 

The results obtained through second approximation 
of the present theory agree well with the exact ones for 
the square plate, at least when E* = 1 (isotropic case). 
It is to be expected, therefore, that the results for 
orthotropic plates, presented in this thesis, will be 
fa,irly accurate. All the results for plates with finite 
length-width ratio indicate that the effects of change in 
the buclfiing pattern must be considered. 

The .'arbitrary parameter a was taken equal to 
square root of (P-Pq)/Pq. This could just as well have 
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been related to the shortening or to the centre deflection. 
For other problems it may be convenient to relate the 
arbitrary parameter to some other property. For example, 
in a thermal buckling problem (Stein, 1959), it is related 
to the ^average rise in temperature beyond that required 
for buckling. 

For experimental verification, to obtain simply 
supported loaded edges, as considered in the present 
investigation, may be impracticable in laboratory. On 
the other hand, it will be easier to subject the panels 
to 'flat-end' loading which results in almost complete 
clamping of the loaded edges. However, if the panel 
tested is long compared with its width, the size and shape 
of the buckles near the centre are expected to be almost 
unaffected by this clamping. Stein (1959) has obtained 
good correlation between theory and experiment for an 
aluminium alloy flat plate having length-width ratio of 
4. Ho experimental results for stiffened plates are 
available . 

Suggestions for Further Work: 

The following problems can be suggested for 


further work: 
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Postbuckling analysis of plates havif^S initial 
geometrical imperfections and/or ini'fci^^ eccentricity 
in loads. 


2 . 


Post buckling strength of plates in sne^r. 
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